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TECHNICAL MEMORANDUM 


THE MSFC/J70 ORBITAL ATMOSPHERE MODEL AND THE DATA BASES FOR 
THE MSFC SOLAR ACTIVITY PREDICTION TECHNIQUE 


1.0 PURPOSE AND SCOPE 

The purpose and scope of this report is to document the MSFC/J70 Orbital 
Atmosphere Model (hereafter referred to as the MSFC/J70 model) and the 13-month 
smoothed values of the 10.7 cm solar radio flux, F 1Q 7 , and the geomagnetic index, 

A , used in the current solar activity prediction technique. 

P 

2.0 GENERAL 

2.1 MSFC/J70 Orbital Atmosphere Model 

The standard, neutral atmospheric density model used at NASA/MSFC for con- 
trol and lifetime studies involving all orbital spacecraft projects, is the MSFC/J70 
Model [1,2]. It is a semi-empirical model that requires values of the 10.7 cm solar 
radio flux and geomagnetic activity index as inputs. 

The model is based on a static diffusion method which defines temperature and 
chemical composition and provides densities in agreement with satellite drag observa- 
tions and, to a somewhat lesser degree, with rocket probe measurements from 90 to 
2500 km altitude. Densities are derived from the empirically determined temperature 
profile and an assumed constant boundary condition at 90 km. Mixing is assumed to 
prevail to an altitude of 105 km and any change in the mean molecular mass below 
this level is assumed to result only from dissociation of oxygen. 

The distribution of mean molecular mass between 90 and 105 km is determined 
empirically in such a way that it results in a ratio of atomic oxygen to molecular 
oxygen of 1.5 at 120 km. All of the recognized variations in the upper atmosphere 
which will be described in the next section, except the ones associated with tidal 
and gravity waves, are included in the model. 

A plot of all the MSFC/J70 constituent concentrations is given in Figure 1, 
versus altitude, for minimum and maximum solar conditions. 

A computational procedure /program listing for the MSFC/J70 Orbital Atmosphere 
Model prediction method is given in Appendix A. The output of the computer pro- 
gram gives temperature, mass density, number density of each constituent, and mean 
molecular weight as a function of altitude (90 to 2500 km), time and location. 

The MSFC/J70 Model evolved from the Smithsonian's Jacchia 1970 (3) model, along 
with two additions from the Jacchia 1971 (4) model. The seasonal-latitudinal varia- 
tions in density below 170 km, and seasonal-latitudinal variations of helium above 
500 km have been incorporated in the MSFC/J70 Model using equations developed by 
Jacchia for his 1971 Thermospheric model. The MSFC/J70 Model is also an integral 
part of the NASA/MSFC Global Reference Atmosphere Model (GRAM) as given in 
Reference 5. 



2.2 Solar Activity Prediction Technique 


Since the MSFC/J70 Model requires values of the 10.7 cm solar radio flux and 
the geomagnetic index as inputs, as well as the orbital parameters, a technique has 
been developed for predicting these values for the time period being considered in the 
various development efforts. This prediction technique is described in detail in NASA 
TM-82462 [6], entitled: "Lagrangian Least-Squares Prediction of Solar Flux (F-.„ „),’’ 

dated January 1984 by R. L. Holland and W. W. Vaughan. 


3.0 DISCUSSION 

3. 1 Atmosphere Model 


Analyses have shown that the following factors contribute to the observed vari- 
ations in the thermospheric (atmosphere above approximately 80 km altitude) density, 
temperature and composition: 

1) Solar activity 

2) Rotation of Earth (diurnal) 

3) Geomagnetic activity 

4) Earth’s rotation about Sun (semiannual) 

5) Seasonal-latitudinal variations of the low thermosphere density 

6) Seasonal-latitudinal variations of helium at all altitudes above 500 1cm 

7) Rapid density fluctuations probably connected with tidal and gravity waves. 

The first six have some regularity and can be modeled with relatively well understood 
degrees of accuracy. Each of these six variations have been taken into account in 
the MSFC/J70 Model. Studies are now underway to incorporate into the MSFC/J70 
Model, representation of the rapid density fluctuations. 

Since temperature plays a minor role in comparison with density, current orbital 
atmosphere models have been developed to represent, insofar as practical and possible, 
the variability of the ambient mass density rather than the temperature. The models 
are based on temperature profiles which have been adjusted so as to produce the 
density values derived from the analyses of satellite orbital decay data. 


3.1.1 Variations with Solar Activity 

In the thermosphere, or orbital atmosphere altitudes, the density is strongly 
influenced by the changing levels of solar activity. The resulting density variations, 
in turn, affect the frictional air drag on orbiting spacecraft. Studies of satellite 
drag data have verified a pronounced correlation between solar activity and density 
variations. However, the reaction of the atmosphere to variations of the Sun is not 
instantaneous. The atmospheric response time is on the order of one day. 
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The ultraviolet solar radiation that heats and causes compositional changes in 
the Earth's upper atmosphere consists of two components, one related to active 
regions on the solar disc and the other to the disc itself. The active-region com- 
ponent varies from day-to-day while the disc component varies more slowly , presum- 
ably with the longer periodicities in the solar activity; i.e, , the approximately 11- 
year solar cycle. The atmosphere has been observed to react in a different manner 
to each of these two components. Jacchia and Slowey (private correspondence) have 
found that the disc component of the solar radiation is, for all practical purposes, 
linearly related to the 10.7 cm solar flux smoothed over six solar rotations (162 days). 
When the short-period oscillations which are caused by the diurnal variation are 
removed, there is essentially an approximate 11-year variation in density that parallels 
the smoothed 10.7 cm solar flux data. 

The night time minimum global exospheric temperature resulting from the varia- 
tion with solar activity can be computed from the following equation [for a geomag- 
netic index (A^) of zero] . 

T C (° K) = 383 + 3.32 F + 1.8 (F - F) (1) 


where: 

_ o 

F = centered 162 day mean F^ ^ cm solar flux (watts/m /cycle /second) 

2 

F = daily mean F^q ^ cm solar flux (watts /m /cycle /second) 


3.1.2 The Diurnal Variation 


Analyses of satellite orbital decay histories have shown that orbital altitude 
atmospheric densities reach a maximum around 2 p.m. local solar time at a latitude 
approximately equal to that of the subsolar point while the minimum occurs between 
3 and 4 a.m. at about the same latitude in the opposite hemisphere. The effect is 
believed to be caused by the absorption of EUV radiation and by heat conduction of 
the neutral gas. The energy that is conducted downward into the lower thermosphere 
and upper mesosphere (altitude range of approximately 50 to 80 km) is mostly lost by 
radiation processes. At an altitude of about 120 km, the time constant for heat loss 
by conduction is on the order of one day or greater. Thus, the diurnal variation is 
not a predominant phenomenon at lower altitudes. Consistency between temperature 
and density cannot always be achieved on a diurnal basis in a quasi-static model; 
therefore, the temperature profile has been used as a parameter so that observed 
density values can be reproduced by the model. Even though the global temperature 
distribution is an artifact developed solely for use in the density and composition 
models, some experimental results are in good agreement. Thom son- scatter tempera- 
ture measurements [7] generally show that the temperature maximum occurs between 
3 and 5 p.m. rather than near 2 p.m. This controversy has not been resolved at 
this time; however, it appears as if there is a phase lag between the density maximum 
and the temperature maximum which cannot be included in the current atmospheric 
models. 

The ratio of day-to-night exospheric (altitude region from 500 to 1000 km) 
temperature was also found not to depend completely upon solar activity although it 
tends to be higher during maximum solar activity and lower during minimum solar 
activity. This ratio has a fairly good correlation with the yearly running mean of the 
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geomagnetic activity index A , but is independent of latitude * The numerical value 
of the ratio varies between 1.27 and 1.40 and has an average of 1.31. 

The correction to be applied to the night time minimum temperature to account 
for the diurnal variation can be computed from the following equation: 


T l (°K) = T c (1+ R sin 2,5 6) (1+ A cos 3,0 (2) 

where : 

T c = night time minimum exospheric temperature (°K) from equation (1) 

R = 0.31 

( 2.5 . 2. 5„\ 

0 - sm 1 \ 

1+ R sin^' & 0 / 

n = \ |lat - ds | 

e = | |lat + ds | 

TAU = HRA -37+6 sin (HRA + 43) 

HRA = hour angle, deg 

LAT = latitude, deg 

DS = declination of Sun, deg 


3.1.3 Variations with Geomagnetic Activity 


Geomagnetic storms usually occur when clouds of charged particles collide with 
the Earth's magnetosphere. These charged particles are believed to be ejected from 
the Sun during the course of a solar flare which is generally a short-lived phenome- 
non. As a result, a large amount of solar radiation is emitted by the flare region 
which subsequently heats the Earth's atmosphere. The heating mechanism is not 
well- understood . 

Analyses of orbital decay histories can give only a blurred picture of the com- 
plex reaction of the upper atmosphere to geomagnetic disturbances. It has been 
shown that the upper atmosphere first reacts in the northern and southern latitude 
polar auroral zones with the energy subsequently propagating toward the equator 
apparently in the form of wave-like perturbations. It appears as if the atmosphere 
reacts with a zero time delay in the auroral zones with the geomagnetic storm effects 
showing up in the equatorial zone about 6 to 9 hr later. An average time lag is 
presently taken as 6.7 hr in the MSFC/J70 Model. 

It is very difficult to adequately include this effect in a quasi- static model. 
The current model calculates the density variations on the basis of a global increase 
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in exospheric temperature. Observations have shown that this is not the case; how- 
ever, for satellite lifetime prediction calculations, an assumed global temperature 
increase is acceptable. For some calculations, such as control dynamics analyses, 
control moment gyro analyses and aerodynamic torques, instantaneous temperature 
(and therefore density), increases for short time periods, and specific locations may 
be required. These studies may require special applications and interpretations of 
the model. Current research is directed toward improving this aspect of the model. 

The geomagnetic activity temperature correction can be computed as follows: 


T„(°K) = 1.0 A + 100 [1 - exp (-0.08 A J] 
G P P 


(3) 


where 

Ap = geomagnetic index value measured 6.7 hr before computation time . 
3.1.4 The Semiannual Variation 


No satisfactory explanation has been found for this variation. It was initially 
assumed that this density variation could be linked with a temperature variation; 
however, data from analyses of more recently orbited satellites showed that the 
original assumption was in error. The amplitude of this density variation is strongly 
height-dependent and variable from year-to-year, with a primary minimum in July 
and principal maximum in October, and a secondary minimum in January, followed by 
a secondary maximum in April. It does not appear to be related to solar activity. 

It has also been found that the semiannual effect varies considerably from solar cycle 
to solar cycle, both in magnitude and in altitude dependence. According to Jacchia 
[8] , the amplitude of this variation is quite large at sunspot maximum but decreases 
toward sunspot minimum. 

The correction to be applied to the temperature to account for this semiannual 
variation can be computed as follows: 


T g (° K) = 2.41+ F [0.349 + 0.206 sin (360 t+ 226.5)] sin (720 t+ 247.6) (4) 


where 




y = 365.2422 days 


DD = day number after January 1 


3.1.5 Seasonal- Latitudinal Variations of the Lower Thermosphere Density 

Presently accepted models assume constant temperature and density values at 
90 km altitude to prevent the models from becoming too complex even though large 
temperature and somewhat smaller density variations are known to exist there. 
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Current models therefore are constructed with a seasonal-latitudinal density variation 
which varies in the vertical from minimal at 90 km, to a maximum at ^ 110 km, and 
then decreases with altitude to 170 km where no appreciable seasonal-latitudinal varia- 
tion has been observed. In the horizontal, the maximum occurs on December 27 at 
the North Pole. These variations are included as additions to the mass densities cal- 
culated as a function of the exospheric temperature used in the model. These varia- 
tions are small and because they occur solely below 170 km altitude they will have 
little effect on orbital lifetime prediction or control requirement calculations. 

After exospheric temperature is computed from 


T e (°K) = Tj^ + T Q + Tg 


(5) 


both local temperature (T^) and mean molecular mass (EM) at altitude, are calculated. 
Finally mass density can be calculated from the following equation: 


DENS (g/cm 3 ) = DENO ex P / 

Z 1 

where 

_ 9 3 

DENO = assumed density at 90 km altitude (3.46 x 10 gm/cm ) 
TO = assumed temperature at 90 km altitude (183PK) 


(EM) G 
(FK) TZ 


dz 


( 6 ) 


MO = mean molecular mass at 90 km altitude (28.878) 
G = gravity (980.665 x [1+ Z/RE] -2 ) cm/sec 2 
Z = altitude, km 

RE = Earth radius (6.356766 x 10 3 km) 

FK = Universal Gas Constant (8.31432 Joule /mole-° K ) 
TZ = local temperature at altitude Z, °K . 


Now the seasonal-latitudinal correction for density can be computed from the following 
equation : 


DENLG(I) = 0.014 (Z-90) exp [-0.0013 (Z-90) 2 ] 


LAT 
LAT | 


sin (2ir$ + 


1.72) x sin 2 (LAT) 


( 7 ) 


where 
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LAT = geographic latitude, deg 
Z = altitude (km) 



$ = semi-annual variation phase ( ft - 36204] / 365 . 2422) 
t = Modified Julian Day 


3.1.6 Seasonal- Latitudinal Variations of Helium 


Experimental results have shown a strong increase r of helium above the winter 
pole. Through an analysis of the intensity of the 10830 A emission line, Bitterberg 
et al. [9] deduced that helium varied seasonally in one year by a factor of 3 to 4 
above 500 km. They pointed out that clear maxima in the helium concentration were 
observed repeatedly in the month of December and January and then were followed by 
a rather steep decline. Minima in the helium concentration were usually observed in 
April and May, but some early ones were observed in March, 

The formation of this helium bulge over the winter pole has been explained by 
a seasonal subsidence of the level at which the diffusion of helium begins. It has 
been shown that a change of this level by 5 km could change the amount of helium 
in the thermosphere by a factor of 2. However, the mechanism of this winter helium 
bulge and its latitudinal dependence are still under investigation . Although the 
mechanism for this migration is unclear ; empirical equations which describe the 
phenomenon are included in the model. These variations influence the computed 
densities only at heights above approximately 500 km. 

The seasonal-latitudinal variation of helium correction can be computed from : 


DLHe = 0.65 


DS 

23.44 


sin 



LAT DS \ 
2 * IDS]/ 


sin 


3 TT 


where 


DS = declination of Sun, deg 
LAT = geographic latitude, deg 


( 8 ) 


3.1.7 Density Waves 


Ambient density waves have been detected throughout the upper atmosphere in 
the height range from 120 to at least 510 km. These fluctuations are believed to be 
caused by tidal and gravity waves. In addition, traveling ionospheric disturbances 
(TIDs) have long been thought of as manifestations of internal gravity waves. Their 
vertical wavelengths apparently increase with altitude until they break down due to 
physical limitations. Density increases on the order of 100 percent have been observed 
to occur over short distances (^10’s of kms) . The waves apparently propagate from 
either south to north, or north to south with maximum horizontal wavelengths on the 
order of 500 to 700 km. Although the current MSFC/J70 model and other models do 
not include variations associated with internally propagating waves, users should be 
cautioned that they have been observed. 
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3.1.8 Total Density Computation 


Total neutral mass density (g/emj is then computed as: 


DL = DLHe + DENLG 


( 9 ) 


where 


DLHe = the log density term modified by Helium S-L variation if greater than 
500 km altitude 

DENLG = the log density S-L component added if less than 170 km alt , 

3 

Finally, the absolute value of density (g/cm ) can be given as: 

DENS - 10 DL . (10) 


3.1.9 MSFC/J70 Model Program Listing and Test Case Example 

The MSFC/J70 Model Fortran program is listed in its entirety in Table 1. It 
represents the current MSFC program running on the HP- 1000. Table 2 contains a 
test case example of total density (kg/m 3) which can be used to compare outputs. 
This test case example was selected from the Space Station Natural Environment 
Design Criteria document, NASA TM- 82585 [10] . It corresponds to the Reference 7 
Table 1 guidance and control density design global mean value of 0.2522 x 10“ H 
kg/m3 at 600 km altitude. Table 3 contains mean orbital total density design values 
for many altitudes from 400 to 1100 km, which can also be used to compare outputs. 
The required inputs for these test cases (Tables 2 and 3) include: 


Date: March 21, 1970 

Time: 1400 UT 

F 10 and F10 = 230 

A^ = 400 
P 

Latitudes = -30° to +30° 

Longitudes = 0° to 350° 

Latitude and longitude increment = 10° . 


3. 2 Solar Activity Data 

The 13-month smoothed values of the 10.7 cm solar radio flux (F lfl and the 
geomagnetic index, A , used as model inputs are listed in Table 4, where they are 

r 

listed as they would be used in a maximum to maximum solar activity prediction. 
Figures 2 and 3 show the 10.7 cm flux and A values, respectively. 

r 
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4.0 CONCLUSIONS 


The MSFC/J70 Orbital atmosphere model provides continuous density, tempera- 
ture, and composition data from 90 to 2500 km for analyses requiring these data. 
Copies of the computer program for this model are available upon request to Chief, 
Atmospheric Sciences Division, NASA, Marshall Space Flight Center, Huntsville, 
Alabama 35812. 

Personnel performing analyses requiring knowledge of any small-scale, short 
time period fluctuations in atmospheric parameters such as those perturbations asso- 
ciated with internally propagating waves should contact the Chief, Atmospheric 
Sciences Division, NASA Marshall Space Flight Center, Huntsville, Alabama 35812, for 
the most current recommended inputs based on results of analyses and studies in 
progress at this time. 
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TABLE 1. MSFC/J70 FORTRAN PROGRAM LISTING (HP-1000 VERSION) 


3; J70MM T” -00 0.04 IS ON CROC 032 USING 0 0151 BLKS R»0000 


0001 

FTN4 

0002 

£EMA 

0003 


0004 

PM :*.+ 

0005 

C+ + 

0006 

C++ 

0007 

C++ 

0008 

C++ 

00 09 

c+ + 

GO I 0 

C+ + 

00 1 1 

C+ + 

0012 

C+ + 

0 0 1 3 

C+ + 

0014 

C+ + 

0015 

C+ + 

O 0 1 6 

C++ 

00 1 7 

c+ + 

00 ! 8 

C+ + 

0019 

c + + 

0 02 0 

c+ + 

0 021 

c+ + 

U 0 i- ci 

C+ + 

0023 

C :■* ■ +■ •+ 

0024 

c 

0025 

C+ + 

0027 


0 028 


0029 


0030 


0031 


0 032 


0033 


0034 

c 

0 035 

C++ 1 

U 036 

c 

0037 


0 038 


0039 

c 

004 0 

f:.+:++. 

0041 

C++ 1 

0042 

c+ ++■ 

.0.043 

c 

0044 


0045 


0046 

c 

0Q4 7 

C++ 1 

0048 

c 

0049 

77*7 

0050 


0051 

1 07 

0052 


0053 


0 054 

1 08 

0055 


0056 

i*’ 

0057 

C 

0 058 

C + + f 


PROGRAM J70MM<3 >>ACI-1 01.281 JACCHIA '70 MODEL < REEDA SYSTEM > 

fc sfc:* * Ht * * * M< + ++ Jfc + +- M + + + :* * * +. * Hc»|c Htrftak & +?+ j* :*:*! :*:* +:+: +.: 4 < 

** 

ESCRlPTiONr PROGRAM ‘ J70MM ' IS THE JAG CHI h '70 MODEL +* 
WHICH WAS CONVERTED FROM THE UN I VAC 1108 ** 

TO THE REEDA SYSTEM. ** 

+.+ 

■4 PUTS .i USER PROVIDES INPUT PARAMETERS VIA 'CRT'. ** 

+•+ 

JTPUTSi A MATRIX GRID OUTPUT IS GENERATED TO THE ** 

HP-2608 PRINTER. A MEAN VALUE OF ALL GRID ** 
VALUES IS ALSO GENERATED, ** 

++ 

4TEST MOD I FICAT I ON i +* 

+■+ 


12 JUL 1985 < CSC > 

.JOHN S. HICKEY <ACI> 533-7590 
M IKE D I CKERSOW < AC I > 533-759 0 

BILL JEFFRIES <CSC> 830-1000 EXT 507 


COMMON /EDATA/ XYR, IDA , MH , IHR, MIN, XM JD, F f 0,F1 0B,GI, ILAT, ILNG, 

, XLAT,NLNG, II 

COMMOW/XDAT A/ AM AT < 73, 37 >,;DENMAT< 73,37 >, DEHL.OG< 73,37 >, HEMAK 7 
HMAT < 73, 37 >,02MAT< 73,37 >, GMAT\ 73 , 37 >, TEMP< 73, 
XN2< 73 , 37 > , WTMAT< 73 , 37 > , EXTEM< 73 , 37 > 

COMMON 7FDATA/ A< 3 00 , 6 > , DE- NLG< 3 0 0 > , D.ENS( 3 0 0 > , DL< 3 0 0 > , EM< 30 0 > , 

TZZ< 300 > , XLATT< 300, XLONG< 3 00 > , Z< 3 0 0 > , XT.EMP< .300 > 


DIMENSION KBUF< 15>, JVARY< 1 1 >, JV< 1 t >, ZALTS< 1 0 > 
DIMENSION IDCBC 276 >, IBUF< 4 0 >, NAME< 3 > , I PAR< 5 > 


CALL RMPAR< I PAR > 
LU m 1 PAR< 1 > 


CALL FTIME< KBUF ') 

WRITE< 6, 1 07 > 

FORMAT< » 1 " > 

DO 7 k'»1 , 1 0 

WR I TE< 6 , 1 08 > < KBUF< I >,1=1,1 5 > 

FORMAT*: ” ********** PROG RAM J7 OMM EXECUTED AT; n ,15A2, 
. ** NhS A/MS FC REEDA SYSTEM *****++*+*” > 

CONTINUE 


U 6 i 



TABLE 1. (Continued) 


0059 

0060 
0 061 
0 062 

0063 

0064 

0065 

0066 
0 067 
0068 

0069 

0070 

0071 
0 072 

0 073 

0074 

0075 

0076 

0077 

0078 

0079 

0080 
0081 
0082 
0033 
0084 
0035 
0086 
0087 
0033 

0089 

0090 

0091 
0032 

0093 

0094 

0095 

0096 

0097 

0098 

0099 

01 00 
0101 
01 02 
01 03 
01 04 
01 05 
0106 

0107 

0108 
01 09 
01 1 0 
01 1 1 
01 12 
01 13 
01 14 
01 15 
01 16 
01 17 
01 IS 


OfcY 

WRITE* LU , 398 > 

398 FORMAT* > 

Cfc2 

WRITE* LU, 399 ) 

3 99 F ORM AT * ” ***** *** ** ** ** * ** ** * * * * * * :+ c :+ t :- fs **** *.* ** ** ** ** ** * * * * * * ** ** ** * 1 / 


»** DESCRIPTION; PROGRAM ' J7GMM ' IS THE JACCHIA '70 **"/ 
"** MODEL WHICH WAS CONVERTED FROM THE **'V 
"** UN I VAC 1108 TO THE REED A SYSTEM, **'V 
"** **"/ 
"** INPUTS! USER PROVIDES ALL INPUT PARAMETERS **”/ 
”** VIA 'CRT-TERMINAL'. **V 
n ** ** 11 / 
”** OUTPUTS! A MATRIX GRID OUTPUT WITH MEAN VALUE **V 
"** IS GENERATED TO THE HP-2608 PRINTER, **V 
M ** **”/ 
"** WRITTEN BY i JOHN $. HICKEY *ACI> 533-7590 **'V 
•'** ** ( V 


. «** MODIFIED BY! BILL JEFFRIES *CSC> 830-1000 EXT 50?**"/ 

, 11 ****** ** ** ****** ************ **** ** ** ****************** *.* 1 ! / 

AO 

C 

Q:h * * ** ** ************ ******* * ** * * ** * * * * * * ** ** **** **** :+r>f« ** 

C** ASK FOR DATE & TIME; YEAR, MONTH, DAY, HOUR, MINUTES ** 

Q*.** ****************************** ** ** .4* *** ** * * **** 

C 

WRITE* LU, 400 > 

400 FORMAT* “^dBENTER^dO Date & Time of Data? * yy ,mm, dd, hh, mm > ; _» > 
READ* LU , * > I YR , MN , I DA , I HR , M I N 

I YR * IYR + 1900 
C 

0**.* ******************** * **** * Hi * * * * * *** 

C+* ask FOR RANGE OF LATITUDE? *-90,90.)** 

C * ***************************** * *** * ***** 

c 

WRITE< LU, 401 > 

401 FQRMAT< 4, ^dBENTERL&dO Range of Latitude? <-90,90; > 

READ< LU , * > ILAT 1 , ILAT2 

C 

C * ** •** ******************** ** * * ***** * * ** * * 

C** ASK FOR LATITUDE INCREMENT?< 5 - 1 0 >** 

C** ************** * * * ********* * *** ****** * * 

c 

WR I TE< LU , 4 02 > 

402 FORMAT< “^dBENTER^xdO Latitude Increment? <5 or 10>; > 

RE AD< LU , * > LINC 

C 

C** COMPUTE NUMBER OF LATITUDE POINTS 
C 

ILAT * < < ILAT2 - ILAT* .VLINO+1 
C 

C*+ STORE LATITUDES INTO XL ATT ARRAY 
C 

XLATT< 1 > * ILAT1 
DO 450 1=2, ILAT 
XLATT< I > « XLATT< I- 1 > + LINC 
450 CONTINUE 

C 

Q* ****************************** *********** 

C++ ASK FOR RANGE OF LONGITUDE ? < 0, 360 > ++ 



TABLE 1. (Continued) 


OF POOR 


quality 


0 1 1 9 
0120 
01 21 
0 1 22 
01 23 
0 1 24 
0125 
0 1 26 
012 ? 
0123 

0129 

0130 
01 3 1 
0132 
01 33 
0134 
01 35 
0 1 36 

0138 
0 1 39 

0140 

0141 
01 42 
01 43 
0 1 44 
0145 
01 46 
0147 
0 1 43 
0143 

0150 

0151 

0152 
01 53 
0154 
01 55 
Cn 56 
01 57 
0 1 58 
0 1 59 
0160 
0161 
0 1 62 

0 1 63 

0164 

0165 
0 1 66 
0167 
0 1 63 
0 1 69 
0 1 70 

0171 

0172 
01 73 
0174 
0 1 75 
01 76 



C‘t -+ + **** ****!* HrHcHeHe * * ** :** ++ + + * * +•+ **.******* 

G 

WRITE<LU,403> 

403 FORMATC "t&dBENTERt&d® Range of Longitude? < .0/36.0? i 
READ< LU , * > I LON 1 , f L0N2 

c 

f: :r tf- +::+■.!+' jfcjf: *•*+! +: * * +< .+ :+ * ++ -Y + :fc :+ + + + +« 

C++ ask for LONGITUDE INCREMENT? <5 - 10>*«* 

Q>t * :+: +:***:+: •+; + +'+ :*»**:♦»:♦« >»«>h 

c 

WRITE< LU, 404 > 

4 0*3 FORMAT< "^dBENTER^^dw Longitude Increment? <5 or 10): 
REAO<LU, + > LLINC 

c 

C** COMPUTE THE NUMBER OF LONGITUDE POINTS 

0 

ILNG « << IL0N2 - I LON 1 >/LL I NO >+ 1 
C 

C++ STORE LONGITUDES INTO XLONG ARRAY 

C 

XLONG< 1 > * I LON 1 
DO 451 1-2/ ILNG 
XLONG< I > = XLONG<I-.1> + LLINC 
451 CONTINUE 

C 

■£:* ^Kit**:**********:**-*:***************************''' ' 

C++ ASK FOR USER INPUTS; 

C++ 2 — ALTITUDE 

C++ FIG - — SOLAR RADIO NO I SB FLUX 

C+ +. F 1 06 — 8 1 -DAY AVERAGE F \ 0 

C++ GI — GEOMAGNETIC ACTIVITY INDEX 

Q . t * ;+i ;+:; * ;4i ********** ** * * * * * * * * * * * * * * ** * * * * Hi * * * * **** 


+:* 

*.+ 

+ * 


6 05 


4 05 
6 06 


WR I TE< LU / 6 05 > 

FORMATS “^dBENTER^^dO Number of’ Altitude Cases? < 1 — 1 0 > s 
READ< LU, + !> MALTS 
DO 606 KK*?/ HALTS 
WR I TE< LU / 4 05 > KK 

FORMAT< "t&dBENTER^Scd© Alti tude< km > # ",12, 

1 11 <2 = 90 to 250 Okm >s > 

RE AD< LU , *0 ZAl.TSC KK > 

CONTINUE 


Cir***#'*#****#*##************************ 

C+ + ASK USER FOR 'INDEX ' VARIABLES ** 
O + II — GEOMAGNETIC INDEX ** 

C++ 12 — DIURNAL EQUATION INDEX +* 

C ?++*:**.*:**************+***** * * * +* + * * * ****** 




« > 


_J* > 


WR I TE< LU , 409 > 

4 09 -FORM AT< dBENTE R dO Geomagnetic Index? <1-KP, 2-AP>* — u > 

READ<LU,*> II 
WRITECLU, 410? 

410 FORM AT< ‘‘^dBENTER^dW Diurnal Equ Index? < 1 -KP , 2-F 1 0B, 3-AVG ? J_ ? 

READ<LU,*> 12 
WRITE<LU,406> 

406 FQRMAT< “£&d8ENTER%&dG» Solar Radio Noise Flux? <F10- 0-460 > ; > 

READ< LU , + > F 1 0 
WR I TE< LU ,407 > 

4 07 FORM AT< ” t&dBENTER^Sc d<? 3 1 -Day Average F10? <F1PB~ .0-250 >; _ > 



TABLE 1 . ( Continued) 


0179 

0180 
0181 
0182 

0 1 33 

0184 

0185 

0186 

0187 

0188 

0189 

0190 

0191 

0192 

0193 

0194 

0195 

0196 

0197 

0198 

0199 

0200 
0201 
0202 

0203 

0204 

0205 

0206 
0207 

02 08 
0209 
021 0 
021 1 
0212 

0213 

0214 

0215 

0216 
0217 
02 IS 
02 1 9 
0220 
0221 
0222 

0223 

0224 

0225 

0226 

0227 

0228 

0229 

0230 

0231 

0232 

0233 

0234 

0235 
02^6 

0237 

0238 


RE AD< LU j * > FI OB 

IRI1 , EG , 2 > WRITE* LU, 408 > 

408 FORMAT* “£&dBENTER| : &d<? Geomagnetic Activity Index? <GI=0~4GG> i > 
IFCri.EQ.1) WRITE* LU, 41 2 > 

412 FORMAT* *%&d8ENTER^&d® Geomagnetic Activity Index? *GI=G-9>; > 

READ* LU , * > GI 
C 
C 

C**# ifcjfc *>♦< ffr * }<f He****** * * >♦* * * ifc tf« Sh * * :4t:+t :+: + :+t :4t»+t % »H H? H* 


c** 

ASK 

FOR ELEMENTS 

TO 

PROCESS <1 - 12 > 

4*4* 

c** 

1 . 

I TEMP 

7. 

IH 

4*4* 

c** 

2. 

IXN2 

8, 

I EM 

4*4* 

C** 

3 « 

102 

9. 

I DEN 

4*4* 

C** 

4, 

10 

10, 

IDLOG 

4*4* 

c** 

5. 

IA 

1 1 . 

IXTEM 

4*4* 

C** 

6. 

I HE 

12. 

ALL OF THE ABOVE 

4*4* 


Q & 4c 4<4« 4«>fi *HH« J* >fc Hofr* 4*4*4* 4*4* * 4* 4*4«4« 4* 4*4* J+t 4* 4*4* 4* 4* 4« 4*4* 4*4* 4* -4*4* 4* 4* 4*4* ifc 


c 

41 1 


c 

c** 

c 

453 

457 

456 


454 

458 

c 

c** 

c 


1 09 


WRITE* LU, 41 1 > 

FORMAT* <72, ,2, 


**4«4< 

SELECT DESIRED ELEMENTS 

TO 

PROCESS ; 

4*4*" ,2, 

"4*4* 

1 , 

TEMPERATURE <K> 

f f 

HYDROGEN *H> 

**“ 

M 4*4* 

2, 

NITROGEN <N2> 

8 . 

MEAN MOLE WT 

1* ,* 

> r ^ 

”4** 

3, 

OXYGEN < 02. > 

9 , 

DENSITY * KG2M3 > 


”4*4* 

4, 

ATOM OXYGEN *0> 

1 0. 

LOG DEN < GM2CM3 > 

*4*” ,2, 


5. 

ARGON * A > 

1 1 . 

EKO TEMP <K> 

4*4* ?l , 2, 


6. 

HELIUM < HE > 

12. 

ALL THE ABOVE 

**",2, 


w 4*4* 4*4*4* 4*4* 4*4* 4* 4* 4* 4* 4* 4* 4*4* 4*4* 4*4* 4* 4c 4*4: 4*4* * 4* 4*4* 4*4* 4*4* 4*4*4* 4* 4* 4* 4*4* 4*4* 4* 4* 4*4* :4<:fs " y 2., 

M t^dBENTER^dO Se lect i ons? * 1 , 3 , 7 > i " > 


READ* LU , * > < JV< J >, J=t , 1 1 > 
INITIALIZE VARIABLE ARRAY 


DO 453 1=1,11 
JVARY* I> = 0 
CONTINUE 

IF* JV< 1 >.NE, 12 > GO TO 456 

DO 457 1*1,11 

JVARY* I > = 1 

CONTINUE 

GO TO 458 

DO 454 1=1,11 

IF* JV< I > . EQ, 0> GO TO 454 

JVARY* JV* I > > = 1 

CONTINUE 

CONTINUE 

PRINT HEADER AND TITLE INFORMATION 


DO 666 KL=1, HALTS 
Z< 1 >=ZALTS* KL > 

WRITE* 6, 1 09 > IYR, MN, IDA, I HR, MIN, ILAT1 , ILAT2, ILON1 , IL0H2, 


Z* O, FI 0, FI OB, G I , I f , 12 


FORMAT * " 1 , 222 , 


PROGRAM NAME; 


J70MM M ,2, 


" DATE & TIME; 

ii ___ , , 

" LATITUDE RNG : 


"14, IK, 12, IK, 12, IK, 12, IK, 12,2, 
V22, 

M47’/M4, 2, 


" LONGITUDE RNG; 


" 22 , 

" * 14, 


14 , 2 , 


14 



TABLE 1, (Continued) 


ORIGINAL PASS' IS 
OF POOR QUALITY 


0239 

024 0 

024 1 

0242 

0243 

0244 

0245 

0246 

0247 

0248 

0249 

0250 

025 1 

0252 C 

0253 C*+- INITIALISE XMJD = 0.0 & IS = t 

0254 C 

0256 XMJD =0,0 

0256 12=1 

025 7 C 

0258 INITIALIZE J COUNTER = 0 

0259 C 

0260 J=0 

0261 C 

0262 C** LOOP TO PROCESS ' XLAT ' DATA 

0263 C 

0264 DO 202 J II«1 , TLAT 

0265 C 

0266 £*■* LOOP TO PROCESS 'XLGMG ' DATA 

0267 C 

0268 DO 201 II =1 , XLNG 

0269 XLAT=XLATT< 1 1 1 > 

0270 XLNG =XLONG< II > 

0271 C 

0272 C*+ LATITUDE ~ LT, LONGITUDE — LD 

0273 0 

0274 DO 200 1=1,12 

0275 J=J+ 1 

0276 C 

027? £*+ CHECK J COUNTER EXCEEDS 30 ? 

0278 C 

0279 IF< J— 30 >52/52,51 

0280 C 

0281 C** RESET J COUNTER = 0 

0282 C 

0283 51 CONTINUE 

0264 J=0 

0285 C 

0266 C*'* CALL 'TME ' SUBROUTINE 

0287 £ 

0288 52 CONTINUE 

0209 CALL TME< MN, IDA, IVR, IHR, MIN, XMJD, XLAT, XLNG, SDA , SNA, DO, DY > 

029 0 C 

0291 Q** CALL # T I HF ' SUBROUTINE 

0292 C 

0294 C+ +• SET RE=0.31++ 

0296 RE = 0,31 

0297 CALL T I NF< F 1 0,F1 OB, GI , XLAT, SDA, SHA, DY, RE, I 1 , I2,TE> 

0298 T =TE 


— “//, 

ALT I TUBE < KM > s ",F5. 1 ,/, 
— . 

FI Os M ,F5 , J j/. 


FI OB? " , F5 . 1 , / , 

">//, 

GIi ",F5.1,/, 

1 — KP , 2—AP : ", 12,/, 

— — ",//, 

1 “KP , 2—F I 08,3-AVGt ”,I2, / 

- ",//> 


15 



TABLE 1* (Continued) 


0299 XTEMP< I > = TE 

0300 C 

0301 C** CALL ' JAC ' SUBROUTINE 

0302 C 

030 3 CALL v#AC< 2< I >, T, T ZZi I >, A< 1 , 1 >, A< 1 , 2 >, A< 1 , 3 >, A< 1 , 4 >, A< I , 5 >, A< 1 , 6 >, 

0304 1 EM< I > , DENS< I >,DL< I >> 

0305 C 

0306 C** INITIALIZE VARIABLES 

0307 C 

03 OS ZZ"Z<I> 

0 309 DUMMY" 0. 

0310 DEN»0, 

0311 DENLG< I >=*0 » 

0312 DUMMY"OL< I > 

0313 DEN»DL< I > 

0314 YOAY=DD 

0315 IF<ZZ-1?0, >20,20,50 

0316 C 

0317 C** CALL 'SLV' SUBROUTINE 

0318 C 

0319 20 CALL SLV< DUMMY, ZZ, XLAT, YDAV > 

032 0 DENLG< I >«DUMMY 

0321 GO TO 40 

0322 50 IF< ZZ-500 » >40,40,30 

0323 C 

0324 C*+ CALL 'SLVH ' SUBROUTINE 

0325 C 

0326 30 CALL SLVH< DEN , A< 1 , 5 > , XLAT, SDA > 

0327 DL< I >«DEN 

0328 40 CONTINUE 

0329 DL< I >«DL< I >*DENLG< I > 

0330 DENS< I > » 10, **DL< I > 

033 1 XLAT"XLAT *< 57 . 2957795 1 > 

0332 C 

0333 C** COMPUTE 0ENMAT< 1 1 , 1 1 1 > 

0334 C 

0335 DENMAT< II , I II > * DENS< I > * 1000. 

0336 TEMP< I I , 1 1 1 > « TZZ< I > 

0337 XN2< I I , II I > « A< I , 1 > 

0338 02MAT< 1 1 , 1 1 1 > « A< 1 , 2 > 

0339 OMAT< 1 1 , 1 1 1 > « A< 1 , 3 > 

0340 AM AT <11, III > » A<I,4> 

0341 HEMAT< II, III > « A<I,5> 

0342 HMAT< 1 1 , 1 1 1 > » A< 1 , 6 > 

0343 WTMAT< I I , I I I > * EM< I > 

0344 DENLOG< I I , I I I > m DL< I > 

0345 EXTEM< I I , I I I > « XTEMP< I > 

0346 200 CONTINUE 

0347 201 CONTINUE 

0348 202 CONTINUE 

0349 C 

0350 C** PERFORM MATRIX PRINTOUT 

0351 C 

0352 I F< JVARY< 9 > , NE . 0 > CALL MA TPR< 1 > 

0353 IF< JVARYC 1 > . NE . 0 > CALL MATPR< 2 > 

0354 I F< J YARY< 2 > , NE . 0 > CALL MATPR< 3 > 

0355 IF< JVARY< 3 > . NE . 0 > CALL MATPR<4> 

0356 I F< JVARY< 4 > , NE . 0 > CALL MATPR< 5> 

0357 I F< JVARY< 5 > . NE . 0 > CALL MA TPR< 6 > 

0358 IF< JVARY< 6 > . NE . 0 > CALL MATPR< 7 > 


16 



TABLE 1. (Continued) 


POOR 


P&QE IS 

quality 


0359 

0360 
036 1 

0362 

0363 

0364 

0365 

0366 

0 367 
0 3.6.8 

0369 

0370 

0371 

0372 

0373 

0374 

0375 

0376 

0377 
0373 
0379 

0330 

0331 

0382 

0383 
0364 

0385 

0386 

0387 

0388 

0389 

0390 

0391 

0392 

0393 

0394 

0395 

0396 

0397 

0398 

0399 

0400 


I F< JVARYC 7 > « HE « 0 > CALL MATPR<9> 

I F< JVARYC 8 > » NE . 0 > CALL MATPR<9> 

I F< OVARYC 1 0 > . HE . 0 > CALL MATPRC 1 0 > 

IF< JVARYC 1 1 NE , 0 > CALL MATPRC 1 f > 

1 04 CONTINUE 
666 CONTINUE 
C 

C** ASK TO CONTINUE? 

C 

WR I TE< LU , 888 > 

888 FORMATC "Process Another Case? <Y or N>: _" > 

READ<LU,609> I PC 
609 FGRMATC A1 > 

I F< I PC ♦ EQ . 1 H Y > GO TO 777 

C 

C** PROGRAM 'J70MM y COMPLETED 

C 

999 STOP 
END 

*EMACXOATA> 

SUBROUTINE MATPRC IP > 

C* + SUBROUTINE 'MATPR ' PERFORMS THE PRINTOUT OF ** 

C++ THE SPECIFIED .MATRIX . ** 

C : + * + ++ ** * sf< >f*Hk * Jftsfci* -ft * 5ft ■+ ; f<tft :ftifs .+ sft :ft :f«.;fe + ++ sft + Jft Jf« *Jft>f«Jf« * Jf« **:**:*::+::** ++ 

C 

C** COMMON STATEMENTS 

C 

COMMON 7EDATA7 I YR , I DA , MH , I HR , M I N , XM OD , F 1 0 , F 1 OB /Cl, I L AT , I LNG , 

. XL AT, XL MG, II 

CQMMGN7XDAT A7 AM AT< 73 , 37 > , OENNATC 73, 37 > , DENLOGC 73 , 37 > , HEMAT < 73 , 37 > , 
HMAT< 73, 37 >, 02MAT< 73 , 37 >,DMAT< 73,37 >, TEMPC 73, 37 >, 

, XH2C 73 , 37 y, WTMATC 73 , 37 >, EXT EMC 73,37 > 

COMMON 7 FDA T A/ AC 300, 6 >, DENLGC 3 0 0 >, OENSC 30 0 > , DLC 300 >, EMC 300>, 

TZZC 300 >, XLATTC 300 >, XLOHGC 300 >, ZC 300 >, XTEMPC 300 > 

C 

C++ DIMENSION STATEMENTS 

C 

DIMENSION I HE ADC 9, 1 1 >, JMTHC 24 > 

C 

C++ DATA STATEMENTS 

C 


04 0 1 

DATA 

JMTH72H JA , 2HN , 2HFE, 2HB , 

2HMA, 

2HR , 

2HAP, 

2 HR , 

2HMA, 2HY 

04 02 
04 03 


2HJU, 2HN ,2HJU,2HL , 
2HN0,2HV , 2HDE, 2HC 7 

2HAU, 

2HG , 

2HSE, 

2.HP , 

2H0C, 2HT 

04 04 

DATA 

I HEAD72HDE , 2HNS , 2HIT, 2HIE 

, 2HS 

, 2H< K 

, 2HG7 

, 2HM3 

,2H> , 

04 05 


2HTE , 2HMP , 2H . , 2HDE 

, 2HG # 

,2H K 

,2H 

,2H 

,2H , 

04 06 


2HCN,2H2>,2H , 2H 

,2H 

,21 H 

,2H 

,2H 

,2H , 

0407 


2HC0,2H2\,2H ,2H 

,2H 

,2H 

,2H 

,2H 

,2H , 

04 08 


2HC 0 , 2H > ,2H ,2H 

,2H 

,2H 

, 2.H 

,2H 

,2H , 

04 09 


2HC ft , 2H j , 2H ,2H 

,2H 

,2H 

,2H 

,2H 

,2H , 

041 0 


2HCH,2HE>,2H ,2H 

,2H 

,2H 

,2H 

,2H 

,2H , 

041 1 


2HCH,2H> ,2H ,2H 

,2H 

,2H 

,2H 

,2H 

,2H , 

04 1 2 


2HME , 2HAN , 2H M , 2H0L 

,2H W 

, 2HT 

,2H 

,2H 

,2H > 

0413 


2HL0 , 2HG , 2HDE, 2HN 

, 2H< G 

, 2HM7 

, 2HCM 

, 2H3 > 

,2H , 

04 1 4 


2HEX , 2H0S , 2HPH , 2HER 

, 2HIC 

, 2H T 

, 2 HEM 

, 2HPC 

, 2HK V 


0415 C 

0416 C++ COMPUTE NGP 
04 1 7 C 
0418 


NPG * ILAT 7 7 



TABLE 1. (Continued) 


0419 

0420 

0421 

0422 

0423 

0424 

0425 

0426 
042? 

0428 

0429 

0430 

0431 

0432 

0433 

0434 

0435 

0436 
043? 

0438 

0439 

0440 

0441 

0442 

0443 

0444 

0445 

0446 
044? 

0448 

0449 

0450 

0451 

0452 

0453 

0454 

0455 

0456 
045? 

0458 

0459 

0460 

0461 

0462 

0463 

0464 

0465 

0466 

0467 

0468 

0469 

0470 

0471 

0472 

0473 

0474 

0475 

0476 
047? 
0479 


C 

O* LOOP FOR HGP TIMES 
C 

IFCMOCKILAT ,7) . NE, 0.1 NPG m NPG + 1 
DO 300 11*1, NPG 
ISTR * < II 
ISTP * ISTR + 6 

IF< ISTP .GT. I LAT > ISTP = ILAT 
C 

C**WRITE HEADER 
C 

INDX1 » <MH-1 >*2 + 1 
JMOH1 » JMTH< INDXI > 

JM0N2 - JMTH< INDXI + I > 

XF<iiIH ,GT • 10> WRITE< 6,7000 > < IHE&CK l , IP), 1 = 1,9), II, JMQNI , JM0N2 
1 IDA, IYR/XNJO, I HR, MIN, 2< 1 >,F1 0,F1 0B,GI, II 

7000 FORMATS 1 H 1 , 29X , 9A2 , 5 1 X , "PAGE “,12, /V 

11X,*DATEi " , 2A2, 12, " ",I4, " JULIAN) MFS.O,” TIME) V2I2, 

2”2 ALT I TUDE< KM > i “ , F7. 1 , IX, M F1 0 s ",F6.2, M F10B? ",F6.2, 

3” Git ** , F6 . 2 , '* O-KP OR 2-AP); Ml) 

IFCMIN.LT, 10) WRITE< 6,8000 > < IHEAD< I, IP), 1=1 ,9>, II, JMOH1 , JM0N2 

1 IDA, IYR, XM JD, I HR, MIN, Z< 1 >/F1 0, FI 06, GI , M 
6000 FORMATC 1H1 ,29H, 9A2, 51 X, "PAGE M27/ 

IfX, "DATE* M , 2A2, 12, n ”,I4,” JULIAN; ",F9.0, " T I ME t ",I2, "0 

2 ”2 ALT I TUDE< KM > i " , F? . 1 , / , IX, "FI 0.* ",F6.2," FI 0B) ",F6.2, 

3" GI I ,, ,F6.2, M <1-KP OR 2-AP>{ ",11 > 

C 

C** PRINTOUT 'XLATT' ARRAY 
C 

WRITE< 6, 7001 > < XLATT< I >, I = ISTR, ISTP > 

7001 FORMAT< 1 HO, 3 OX, "< -SOUTH) LATITUDES < +N0RTH > V1X, "LOW* M , 8X, 
1?<F5.0,5X>> 

UR2TE< 6/7003 > 

7003 FORMATC 1H , "<■— WEST W1X, "< +EAST)"/) 

G 

C** PRINTOUT 'XLONG ' ARRAY 
C 

DO 250 J»1 , ILNG 

IF< IP ,EQ . 1 ) WR I TE< 6 , 7 0 02 > XL ONG< J > , < DEHM AT< J , I > , I * I STR , I STP > 
IF< IP.EQ . 2 1 WRITE< 6, 7002 > XLONG< J>,< TEMP< J, I >, I*ISTR, ISTP) 

I F< I P , EQ , 3 ) WRITE< 6 , 7002 > XLONG< J),< XN2< J , I ) , I » I STR , ISTP!) 
IF< IP . EQ . 4 > WRITE< 6, 7002 > XLOHG<J>,< 02MATC J, 2 >, I»ISTR, ISTP ) 
IF< IP ,EQ ■ 5) WRITE< 6,7002) XLOHG< J>,< OMAT< J , I >, I-ISTR, ISTP) 
IF< IP , EQ , 6 > WRITE< 6, 7002 > XLONG< J >, < hMATC J, I >, I-ISTR, ISTP ) 
IF< IP.EQ.?) WRITE< 6 , 7002 > XLONG<J>,< HEMATC J, I >, I-ISTR, ISTP ) 
IF< IP . EQ • 8 > WRITEC 6, 7002 ) XLONG< J>,< HMATC J, I >, I»ISTR, ISTP) 
IF< IP.EQ, 9 > WR I TE< 6,70 02 > XLONG<J>,< WTMATC J, I >, I -I.STR, ISTP > 
IF< IP .EQ .10) WRITE< 6, 7002 > XLONG< J >, < DENLOG< J, I >, I*ISTR, ISTP > 
IF< IP.EQ. 1 1 > WRI TE< 6 , 70 02 > XL0NG<J>,< EXTEM< J, I ), I “ISTR, ISTP > 

7002 FORMAT< 1 H ,F4 . 0, 4X, 7EI 0 . 4 > 

250 CONTINUE 

300 CONTINUE 
C 

C** COMPUTE AND PRINT MEAN 
C 

NPTS«ILNG*ILAT 
TOT-O *0 

DO 10 1*1 , ILAT 
DO 20 J»1,ILNG 

IF< IP.EQ, 1 > TOT*DENMAT< J, I >+TOT 



TABLE 1. (Continued) 


nr r.^ ' * B§ 

0F POOR QUALITY 

0479 IRIP.EQ.2) TOT- TEMP< J, I >+T.0T 

0480 IF< IP • EQ . 3 > TOT* XN2< J, I >+TOT 

0481 IF< IP »EQ 4 4 > TOT* 02MAT< J, I >+TOT 

0482 IF< IP . EQ . 5 > TOT* GNAT < J , I >+TGT 

0483 IF< IP. £G.6> TOT* ANAT< U, I >+TGT 

0484 IF< IP , EQ , 7 > TOT* HEMAT< 0, I >+TQT 

0485 IF< IP , EQ . 8 > TOT* HMAT< J,1 >+TGT 

0486 IF< IP * EQ . 9 > TOT* WTMAT < J , I >+TOT 

0487 IF< IP »EQ . 1 0 >TGT*DENLOG< J, I >+TGT 

0488 IF< IP.EQ , 1 1 >T0T* £XTEM< J, I >+TGT 

0489 20 CONTINUE 

0490 1 0 CONTINUE 

049 1 XMEAN*TOT/FLOhT< NPTS > 

0492 WR I TE< 6 , 444 > NPTS , XNEAN 

0493 444 FQRMAT< //, IX, "Humber of Da t-a Valuer i M 5, 5X, " Mean Value* ", 

0494 +E1 0,4 > 

0495 C 

0496 C++ RETURN TO CALLING PROGRAM 

0497 C 

0498 RETURN 

0499 END 

050 0 SUBROUT I NE T I NF< F 1 0 , F 1 0B , G I , XL AT , SDA , SHA , O Y , RE , 1 1 , 1 2 , TE > 

050 1 C+ + + ++++++ + + +++ + + + + •+ + + + + ++ + + + + + + + + ’++++++++++++ ++++ ; ++ : +++*+ ++++++ 

0502 C++ SUBROUTINE 'TINF ' CALCULATES THE EXOSPHERIC TEMPERATURE ++ 


0503 C++ ACCORDING TO JACCHIA SAG NO, 313, 1970, ** 

0504 C++ ++ 

0505 C++ F 1 0 » SOLAR RADIO NOISE FLUX < XE-22 WATT3/M++2 > ** 

0506 C++ FI 0B* 81 -DAY AVERAGE F10 ++ 

0507 C++ GI » GEOMAGNETIC ACTIVITY INDEX ++ 

0508 C++ LAT * GEOGRAPHIC LATITUDE AT PERIGEE < IN RAD > ++ 

0509 C++ SDA * SOLAR DECLINATION ANGLE < IN RAD;) ** 

0510 C++ SHA * SOLAR HOUR ANGLE ++ 

0511 C++ DY * O/Y < DAY HUMBER/ TRGP I CAL YEARX; 1 ++ 

0512 C++ II * GEOMAGNETIC EQUATION INDEX* 1 ~GI*KP, 2— GI*AP. ** 

0513 C++ 12 * DIURNAL EQU INDEX < 1 ~~R< KP >, 2— R< FI OB >, 3-R< AVG >, ** 

0514 C++ RE * DIURNAL FACTOR KP, Ft 08, AVG* ++ 

051 5 C++ ++ 

0516 C++ CONSTANTS — C*SOLAR ACTIVITY VARIATION. •+■+ 

0517 C++ ~ BETA, ETC. * DIURNAL VARIATION. ++ 

0518 C++ — D=GEQMAGNET I C VARIATION. ++ 

0519 C++ E*SEM I A NN UAL VARIATION, ++ 

0520 C++ ++ 


052 1 Q.+++ +++++++++++++•+++++ ++++++++++++++++++++++++ ++++++++++++++++ 

0522 C 

0523 C++ DATA STATEMENTS 

0524 C 

0525 DATA C 1 /383 . 0/ 

0526 DATA C2/3 , 32/ 

0527 DATA C3/1 .80/ 

0528 DATA PI/3.14159/ 

0529 C 

0530 C++ PERFORM CALCULATIONS 

053 1 C 

0532 CON « PI/180. 

0533 BETA* -37.Q+C0N 

0534 GAMMA* 43 . Q+COM 

0535 P * 6 . O+CON 

0536 XM » 2,5 

0537 XNN * 3.0 

0-538 C 
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0539 

0540 

0541 

0542 

0543 

0544 

0545 

0546 

0547 
0543 

0549 

0550 

0551 

0552 

0553 

0554 

0555 

0556 

0557 
0553 

0559 

0560 

0561 

0562 

0563 

0564 

0565 

0566 

0567 

0568 

0569 

0570 

0571 

0572 

0573 

0574 

0575 

0576 

0577 

0578 

0579 

0580 

0581 

0582 

0583 

0584 

0585 

0586 

0587 

0588 

0589 

0590 

0591 

0592 

0593 

0594 

0595 

0596 

0597 

0598 


C+* INITIALIZE GEOMAGNETIC VARIATION VARIABLES 
C 

D1 * 28.0 

D2 - 0.03 

03 * 1.0 

04 - 100.0 

05 * -0.08 

C 

C++ INITIALIZE SEMIANNUAL VARIATION VARIABLES 
C 


El - 

2.41 

£2 « 

0.349 

E3 * 

0.206 

E4 » 

360 . *CON 

E5 » 

226, 5+CON 

E6 * 

720.+CON 

E7 - 

247,6+CGN 

E8 - 

0.1145 

E9 ■« 

0,5 

El 0* 

E4 

Ell* 

342,3+CON 

El 2* 

2/16 


C 

C++ SOLAR ACTIVITY VARIATION 
C 

TC Cl + C2+F10S + C3*<Ft 0 - Ft OB.} 

C 

C++ OIURNAL VARIATION 
C 

ETA - 0 . 5*ABS< XLAT - SOA> 

THETA * 0 * 5+ABS< XLAT + SPA) 

TAU * SHA + BETA + P*SIH<SHA + GAMMA) 
TPI«2+PI 

IF< TAU > 21 0/230,230 
210 IF< TAU+PI > 220,250/250 
220 TAU“TAU+TPI 
GO TO 210 

230 IF< TAU-PI > 250,250,240 
240 TAU«TAU-TPI 
GO TO 230 
250 CONTINUE 
C 

£*** h«+ ++++ + + + + * ++ ++ + + + * + + ++ 

C++ SET R » RE CRE » 0.31 > *♦ 

Qrtofcjf' 

c 

RE « 0.31 
R ** RE 

A1 »< S I N< THE T A > >**XM 
A2 «< COS< ETA > )«XM 
A3 “< COS< TAU72 . > >++XNN 
B1 » 1.0 + R+A1 
B2 *< A2-A 1 V< 1. + R+A1 > 

TV » 81 *< 1. + R+B2+A3) 

TL * TC+TV 
C 

C++ GEOMAGNETIC VARIATION 
C 

IF < It— t > 50,50,60 
50 TG * D1+GI + D2+EXP< G I ) 



TABLE 1. (Continued) 


0599 GO TO 70 

0600 60 TG « 03+GI + D4+< T -EKP< 05+GI > > 

0601 70 CONTINUE 

06 02 C 

0603 C++ SEMIANNUAL VARIATION 

0604 G 

0605 G3 « 0 , 5*< 1 . 0 + SIN<E1 0+DY +£t 1 > > 

0606 G3 - G3++E12 

0607 TAU1 - DY + E8+< G3 - E9 > 

0608 G1 * E2 + E3*<SIN<£4*TAU1 + E5>> 

0609 G2 » S I N< E6-+TAU1 + E7> 

0610 TS « El + F1GB+G1+G2 

0611 C 

0612 C++ EXOSPHERIC TEMPERATURE 

0613 C 

0614 TE = TL + TG + TS 

0615 C 

0616 C++ RETURN TO CALLING PROGRAM 

0617 € 

0618 RETURN 

0619 END 

0620 SUBROUTINE TMECMN, IDA, I YR, IHR, MIN, NM JO, XLAT/XLHG, SOA, SHA, DO, DY > 

062 1 £++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 


0622 C++ SUBROUTINE ' THE ' PERFORMS THE CALCULATIONS OF THE SOLAR DECLI- ++ 

0623 C++ NATION ANGLE AND SOLAR HOUR ANGLE. ++ 

0624 C++ ++ 

0625 C++ INPUTS: MN a MONTH ++ 

0626 C++ IDA = DAY ++ 

0627 C++ IYR = YEAR *+ 

0628 C++ IHR = HOUR ++ 

0629 C++ MIN = MINUTE ++ 

0630 C++ XMODa MEAN JULIAN DATE < IF~0 THEN XMJD IS CALCULATED > ++ 

0631 C++ XLATa LATITUDE < INPUT-GEOCENTRIC LATITUDES ++ 

0632 C++ XLNQa LONGITUDE* INPUT-GEOCENTRIC LONGITUDE, - 1 80, +1 80 > ++ 

0633 C++ ++ 

0634 C++ OUTPUTS * SDA « SOLAR DECL INATION ANGLE < RAD !> ++ 

0635 C++ SHA = SOLAR HOUR ANGLE < RAD > ++ 

0636 C++ DD a DAY NUMBER FROM 1 JAN. +* 

0637 C++ DY * DD/TROPT CAL YEAR 


0638 C+++++ +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 

0639 C 

0640 C++ DIMENSION STATEMENTS 

0641 C 

0642 DIMENSION IDAY< 12 > 

0643 C 

0644 C++ DATA STATEMENTS 

0645 C 

0646 DATA IDAY/31 , 28, 3 1 ,30, 31 , 30, 3 1 , 31 , 30, 31 , 30, 31 / 

0647 DATA YEAR/365,2422/ 

0648 C 

0649 C++ SET CONTSTANTS 

0650 C 

0651 XLATaXLAT/57 .29577951 

0652 YR«IYR 

0653 JaIYR~4+< IYR/4 > 

0654 IF< J>1 0*5, 10 

0655 5 IDAY<2>»29 

0656 10 CONTINUE 

0657 IF < MN— 1 > 3, 3,4 

0653 3 00*= I DA 
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(Continued) 


0659 

0660 
0661 
0662 

0663 

0664 

0665 

0666 
066? 
0668 

0669 

0670 

0671 

0672 

0673 

0674 

0675 

0676 

0677 

0678 

0679 

0680 
0681 
0682 

0683 

0684 

0635 

0636 
0687 
0638 

0689 

0690 

0691 

0692 

0693 

0694 

0695 

0696 

0697 

0698 

0699 

0700 
0 ? 01 

0702 

0703 

0704 

0705 

0706 

0707 

0708 

0709 

0710 
071 1 

0712 

0713 

0714 

0715 

0716 

0717 

0718 


GO TO 6 
4 KE-MN-1 
ID*0 

DO 20 I»1,KE 
ID - ID + IDAYCI) 

20 CONTINUE 

ID * ID + IDA 
0D » ID 
6 DY * DD/YEAR 
C 

C** COMPUTE MEAN JULIAN DATE IF XMJD = 0 
C 

IF<XMJD> 30,25,30 

25 XMJD - 2439856. + 365 . *< YR-1 968 . > + DD 
LDD » CIYR-1965V4 
XMJD * XMJD + LDD 
30 FMJD » XMJD - 2435839. 

C 

C** COMPUTE GREENWICH MEAN TIME IN MINUTES GMT 
C 

XHR * I HR 
XMIN = MIN 
GMT * 60*XHR + XMIN 
C 

C** COMPUTE GREENWICH MEAN POSITION - GP < IN DEG) 

C 

XJ » < XMJD - 2415020. 0V< 36525. 0) 

At ®99, 6909833 
A2 « 36000.76854 
A3 « 0.00038708 
A4 « 0.25068447 

GP « A1 + A2*XJ + A3*XJ*XJ + A4*GMT 
N * GPS360, 

XN * N 

GP » GP - XN*360. 

C 

C** COMPUTE RIGHT ASCENSION POINT - RAP < IN DEG > 

C 

C** 1ST CONVERT GEOCENTRIC LONGITUDE TO DEG LONGITUDE - WEST NEG + EAST 
C 

I FACT » XLNG/180. 

XFACT * I FACT 
XLNG 3 *XLNG-36 0 . *XFACT 
C 

RAP * GP + XLNG 
N * RAP/360. 

XN * N 

RAP m RAP - XN*360. 

C 

C** COMPUTE CELESTIAL LONGITUDE - XLS < IN RAD) ZERO TO 2PI 

C 

B1 * 0.017203 
82 - 0.0335 
B3 - 1.410 
Y1 » B1*FMJD 

XLS * Yt + B2*SIN<Y1> - B3 
TP I - 6.28318 
N » XLS/TPI 
XN » N 

XLS - XLS - XN*TPI 
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0719 

0720 

0721 

0722 

0723 

0724 

0725 

0726 

0727 

0728 

0729 

0 730 

0731 

0732 
0 733 

0734 

0735 

0736 

0737 

0738 

0739 

0740 

074 1 

0742 

0743 

0744 

0745 
07 4.6 

0747 

0748 

0749 

0750 

0751 

0752 

0753 

0754 

0755 

0756 

0757 

0758 

0759 

0760 
076 I 

0762 

0763 

0764 

0765 

0766 

0767 

0768 

0769 

0770 

0771 

0772 

0773 

0774 

0775 

0776 

0777 
0773 


C 

C +* COMPUTE SOLAR DECLINATION ANGLE - SDA < IN RAD > 

C 

B4 = <TPI/360, >*23.45 

SDA » AS I N< S I N< XLS >*S I H< 84 > > 

C 

C++ COMPUTE RIGHT ASCENSION OF SUN - RAS < IN RAD > ZERO TO 2PI 

C 

RAS » AS IN< TAN< SDA .V*TAM< B4 > > 

C 

C++ PUT RAS IN SAME QUADRANT AS XLS 
C 

PI * 3 . 14159 
PI2 * PI/2, 

PI 32* 3.*PI2 
RAS * A8S< RAS > 

TEMP « ABS< XLS > 

IF< TEMP - PI2> 130,130,100 
1 00 IF< TEMP - PI > 105,1 05, 1 1 0 
1 05 RAS » PI - RAS 
GO TO 130 

110 IRTEMP - PI 32 > 115,115,120 
1 15 RAS * PI + RAS 
GO TO 130 

120 RAS * TP I ~ RAS 
130 IF < XLS > 135,140,140 
135 RAS = -RAS 
140 CONTINUE 
C 

C++ COMPUTE SOLAR HOUR ANGLE - SHA < IN OEG > - - 

C 

SHA * RAP-*K PI / 1 80 . > - RAS 
IF<SHA> 210,230,230 
210 IF< SHA+PI > 220,250,250 
220 SHA* SHA+TP I 
GO TO 210 

230 IF< SHA-PI > 250,250,240 
240 SHA* SHA— TP I 
GO TO 230 
250 CONTINUE 

C 

C** RETURN TO CALLING PROGRAM 


C 

IDAY<2>*23 

RETURN 

END 

SUBROUTINE 0AC< Z, T , TZ, AN, A02 , AG, AA , AHE , AH, EM, DENS, DL > 

C ++ Jf^ili**^**^*^^*:* + ■¥ *•- tf< :+c:+t :+c*:fc.4<:+c An sfr % sfc :+::+c * * ** ** sfc* + + 

C** SUBROUTINE ' JAC ' PERFORMS THE SIMPSONS RULE QUADRA- ** 


C++ TURE < SRQ4 > IMPLEMENTED BY G . F. KUNCIR, ** 

C++ ++ 

C++ A * LOWER LIMIT OF INTEGRATION ** 

D * UPPER LIMIT OF INTEGRATION ** 

C++ FUNC * INTEGRAND FUNCTION SUBPROGRAM ** 

C++ EPS « RELATIVE ERROR CONVERGENCE CRITERION ** 

C++ M * MAXIMUM NUMBER OF INTEGRATIONS ** 

C++ R » RESULT OF INTEGRATION ■++ 

C++ N * NUMBER OF INTEGRATIONS REQUIRED TO FIND R ** 

C++ ++ 





TABLE 1, (Continued) 


0779 

0700 

0781 

0702 

0783 

0784 

0785 

0786 

0787 

0788 

0789 

0790 

0791 

0792 

0793 

0794 

0795 

0796 

0797 

0798 

0799 

0800 
0801 
0802 

0803 

0804 

0805 

0806 

0807 

0808 
0809 
081 0 
081 1 
0812 

0813 

0814 

0815 

0816 

081 7 

0818 

0819 

0820 
0821 
0822 

0823 

0824 

0825 

0826 

0827 

0828 

0829 

0830 

0831 

0832 

0833 

0834 

0835 

0836 

0837 

0838 


C 

C++ DIMENSION STATEMENTS 
C 

DIMENSION ALPHA< 6 >, EI< 6 >,DI< 6>/B< 7 >,DIT< 6 > 

C 

C++ DATA STATEMENTS 
C 

DATA QQ7 100.7 

DATA ALPHA/0,0,0. 0, 0,0, 0,0,-. 380, 0.07 
DATA AV76 . 02257E237 

DATA EI/28.0134,31 . 9988, 1 5 . 9994, 39 . 948, 4 . 0026, t , 00797/ 

DATA 8728.15204,-0.085586, 1 . 2840E-04, -1 . 0056E-05,-1 . 021 0E-05, 
1 . 5044E-06, 9 . 9826E-G8/ 

DATA QN7.781107 
DATA QQ2/.209557 
DATA QA7 . 0093437 
DATA QHE71 .289E-057 
DATA FK78.314327 

C++ 

C++ SET VARIABLES 
C++ 

ALPHA< 1 > * 0.0 
ALPHA<2> = 0.0 
ALPHA<3> * 0.0 
ALPHA< 4 > = 0.0 
ALPHA< 5 > -.38 

ALPHA<6> * 0.0 
AV = 6, 02257E23 
.EK'D » 28.0134 
EI< 2 > * 31.9988 
E I < 3 > = 15.9994 
EI< 4> » 39.948 
EI< 5 > » 4.0026 
EI<6> » 1 .00797 
B< 1 > * 28. 15204 
B< 2 > .* — 0 . 085586 
B< 3 > « t . 284QE-04 
B<4 > « -1 . 0056E-05 
B< 5 > « -1 , 0210E-05 
B< 6 > = 1 . 5044E-06 
B< 7 > » 9.9826E-08 
ON « .78110 
Q02 - .20955 
QA « .009343 
QHE - 1 .289E-05 
FK - 8.31432 
C 

C++ PERFORM CALCULATIONS 
C 

TI » T 

TX«444.3807+. 02385+TI-392 , 8292+£XP< - , 0021357+TI > 

A2=2.*< T-TX >73. 14159265 
DIT< 6 >=*0 . 

M«f 0 

EPS 6 * .0001 

T 1 « 1 . 9*< TX- 1 83 . >735 . 

T4=3 . *< TX- 183.-2 . +T 1 +35 . 73 . >7< 35 . ++4 > 

T3=-T1 7< 3 , +35 . ++2 >+4 . +T4+35 . 73 . 

TZ»TX+T1*<Z-125, >+T3*<2-125. >**3+T4*< 2- 1 25 , >**4 
IF <2-105. > 43,43,40 
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0839 

0840 

0841 

0842 

0843 

0844 

0845 
0848 

0847 

0848 

0849 

0850 

0851 

0852 

0853 

0854 

0855 

0856 

0857 

0858 

0859 

0860 
0861 
0862 

0863 

0864 

0865 

0866 

0867 

0868 

0869 

0870 

087 1 

0872 

0873 
■0074 

0875 

0876 

0877 

0878 

0879 

0880 
0881 
0832 

0883 

0884 
0835 
0886 

0887 

0888 

0889 

0890 

089 1 

0892 

0893 

0894 

0895 

0896 

0897 

0898 


43 22 * 2 “ Qw 

EM“B< 1 >+B< 2 >+Z2+B< 3 :>+Z2+*2+B< 4 >*Z2**3+B< 5 >+22++4+B< 6 >+Z2**5 
1 +B< 7 >*Z2**6 
D=2 

70 CONTINUE 
8=90 , 

FA=B< 1 >+B< 2 >*< A-QQ >+B< 3 >+< A-GQ >**2+B< 4 >*< A-QQ }+*3+B< 5 >*< A-QQ >**4 
1 +B< 6 >*< A-QQ >**5 +B< 7 >+< A-QQ }**6 
FA=FA*9 , 80655/< < 1 , +A/6 , 35 6766E+3 >**2 > 

FA=FA/< TX+T1 *< A- 1 25 , >+T3+< A- 125 . >**3 +T4*< A- 1 25 , >**4 > 

FD=B< 1 >+B< 2 >+< O-GQ >+8< 3 >+< D-GQ >**2+B< 4 >*< O-QQ >**3+B< 5 .*>+< D-QG >**4 
1 +B< 6 >*< D-QQ>**5 +B< ? >+< O-QQ >+*6 
FD=FD*9 ,8 0665/< < 1 , +0/6, 356766E+3 >*+2 > 

FO=FD/< TX+T 1 +< D- 125 , >+T3+< 0- 1 25 , >**3 +T4+<D~125. >**4> 

C 

C++ INITIALIZE COUNTERS 
C 

N»0 

PREV=0. 

SONE=< O-A >+< FA+FD V2 , 

71 N=N+1 

IF < N-M > 72,72,75 

72 NINT=2+*N 
STUO=0. 

DEL=< D-A >/FLOAT< HINT > 

DO 73 1 = 1 , NINT, 2 
X«A+DEL*FLOAT< I > 

FX=B< 1 )+B< 2 >*< X-GQ >+B< 3 >+< X-QQ >**2+B< 4 >*< X-GQ >++3+B< 5 .’>+< X-QQ >**4 
1 +B< .6 >*< X-QQ >++5 +B< 7 >+< X-QQ >++6 
FX=»FX+9 . 8 0665/< < 1 . +X/6 . 356766E+3 >**2 > 

FX=FX/<TX+T1+<N-125. >+T3*CX-1.25 . >**3 -+T4*< X-1 2.5, >*+4 > 

73 STWO=STWO+FX 
CUR=SONE+4 , +DEL+STW0 

IF < EPS*ABS< CUR >-ABS< CUR-PREV > > 74,75,75 

74 PREV*CUR 

SONE»< SONE+CUR V4 . 

GO TO 71 

75 R=CUR/3 

IF < Z— 1 05 . > 44,76,44 
44 IF <0-1.05*. > 76,55/76 

76 DENS =3 , 46E-9+ 1 83 , +EM+EXPC -R/FK >/.< TZ+28 . 878 > 

DL=ALOGT< DENS > 

P AR=A V +DENS/EM 
AN=ALO£T< QN+EM+P AR/2S , 96) 

AA=ALOGT< Q A+EM * PAR/ 28 . 96 > 

AHE=ALOGT< QHE+EM+P AR/28 . 96. > 

AO=ALOGT< 2 . *PAR+< 1 , -EM/28 , 96 > > 

A02=AL0GT< PAR*< EM+< 1 . +002 V28 , 96-1 . Y> 

AH=-0. 

C 

C++ RETURN TO CALLING PROGRAM 
C 

RETURN 

C 

C++ CONTINUE CALCULATIONS 
C 

40 23=1 05, 

TZ3=TX+T1 +< 23-125 . ,HT3*< 23- 1 25 >**3+T4+< 23-1 25 >++4 
2M3=B< 1 >+B<2>* 5 . +B< 3 >+ 25.+B<4>+ 125,+B<5>+ 5 , +*4 . +B< 6 >* 5.++S. 
1 +B< 7 >+ S.++6. 
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0899 



D® 1 05 « 

0900 



GO TO 70 

0901 

55 


DENI ®3 . 46E-9*t 83 . *ZM3*EXP< -RZFK >/< TZ3*28 , 878 > 

0902 



PftR»aV* 0EN1/ZM3 

0903 



DI< 1 >»QN*ZM3*PARX28,96 

0904 



DI<2 >«PAR*<ZM3*< 1 , +Q02V28 , 96-1 . > 

0905 



DI< 3 >®2 . *PAR*< 1 , -ZM3/28 . 96 > 

0906 



DI< 4 >®QA*ZM3*PAR/2S , 96 

0907 



DI< 5 >»QHE*ZM3*PAR/28 . 96 

0908 



IF< Z— 125 . > 56,56,90 

0909 

56 


CONTINUE 

0910 



A 1 « 1 05 . 

0911 



FA 1 *9 . 80665/< < 1 . +A 1 /6 . 356766E+3 >**2 > 

0912 



FA1«FA1/<TX+T1*< A1-125, >+T3*< A 1 - 1 25 . }**3+T4*< A 1 - 1 25 . !>**4 > 

0913 



D1»Z 

0914 



FD1 *9 . 30665/< < 1 . +D1 /6 . 3567 66E+3 >**2 :> 

0915 



IF< D1 — 1 25 , > 45,45,50 

0916 


45 

FD 1 ®FD 1 /< TX+T 1 »*< D 1 — 1 25 « >+T3*< D1 -1 25 . >**3+t4*<01 -125 ♦ >**4> 

0917 



GO TO 51 

0918 

50 


FD1»FD1/<TX+A2*ATAN<T1*<D1-125. >*< 1 . +4 . SE-6*< D1 -1 25. >**2 . 5 >/A2 > > 

0919 



TZ»TX+A2*ATAN< T 1 *< Z- 1 25 . >*< 1 , +4 , 5E~6*< Z- 1 25 . >**2 , 5 VA2 > 

0920 

51 


N®0 

0921 



PREV»0 

0922 



SONE«<D1-A1 >*< FA1 +FD1 V2, 

0923 

81 


N«N+1 

0924 



IF <N-M> 82,82,85 

0925 

82 


NINT»2**N 

0926 



STWO®0, 

0927 



DEL»<D1-A1 VFL0AT< NINT > 

0928 



DO 83 I®1 , NINT, 2 

0929 



XI »A1 +DEL*FL0AT< I > 

0930 



FX1»9,80665/<< 1 , +X1 X6 , 356766E+3 >**2 > 

0931 



IF< XI —125 ♦ > 46,46,52 

0932 


46 

FX 1 ®FX 1 /< TX+T 1 *< XI — 1 25 . >+T3*< XI -1 25 . >**.3+T4*< X 1 - 1 25 , >**4 > 

0933 



GO TO 83 

0934 

52 


FX 1 «FX 1 X< TX+A2*AT AN< T 1 *< X 1 — t 25 , >*< 1 . +4 . 5E“6*< X 1 ~ 1 25 . >**2,3VA2» 

0935 

83 


STUK3*STWO+FX1 

0936 



CUR»SONE+4 , *DEL*STUO 

0937 



IF < EPS*ABS< CUR >“ABS< CUR-PREV > > 84,85,85 

0938 

84 


PREV=CUR 

0939 



SONE-< SONE+CUR >X4 , 

0940 



GO TO 81 

0941 

85 


R«CUR/3 , 

0942 



DO 41 1*1,5 

0943 



DIT< I )»D I < I >*< TZ3XT2 >**< 1 . +ALPHA< I > }*EXP< ~EI< I >*RXFK > 

0944 


41 

CONTI NIC 

0945 



D£NS®0 

0946 



DO 42 1*1,6 

0947 



DENS*DEN3+EI< I >*DTT< I VAV 

0948 

42 


CONTINUE 

0949 



EM®OENS*AVX<DXT< 1 >+DIT< 2 >+DIT< 3 >+DIT< 4 >+DIT< 5 >*DIT< 6 >> 

0950 



DL«ALOGT< DENS > 

0951 



AN ®ALOGT<DXT< 1 >) 

0952 



A02*AL0GT<DXT<2>> 

0953 



AO =ALOGT<DIT<3>> 

0954 



AA «ALQGT<DXT<4>> 

0955 



AHE«*AL0GT< DIT< 5 > > 

0956 



IF< Z-50G , > 47,48,48 

0957 


47 

DIT< 6 >»1 0 . *M*< —6 > 

0958 


48 

AH®AL0GT< DIT< 6 >> 



TABLE 1. (Continued) 


0959 

0960 

0961 

0962 

0963 

0964 

0965 

0966 

0967 

0968 

0969 

0970 

0971 

0972 

0973 

0974 

0975 

0976 

0977 

0978 

0979 

0980 

098 1 

0982 

0983 

0984 

0985 

0986 

0987 

0988 

0989 
099 0 

099 1 

0992 

0993 

0994 

0995 

0996 

0997 

0998 

0999 
1 000 
1 001 
1 002 
1 003 
1 004 
t 005 
1 006 
1 007 
) 008 
1 009 
1 01 0 
1011 
1 012 
1 013 
1 0t4 
1 015 
1 016 
1017 
1 01 8 


AN =AMAX1 < -0 . , AN > 

A02«AMAX1<-0. ,A02> 

AO **AMAX1 < ~0 . , AO > 

A A =AMAX 1 < —0 , , A A > 

AHE=AMAX 1 < -0 . , AHE > 
ah =AMAX1 < — 0 . , AH :> 
e 

C++ RETURN TO CALLING PROGRAM 
C 

RETURN 

e 

C++ CONTINUE CALCULATIONS 
C 

90 S=TX+A2*ATAN< T 1 +375 . •+< 1 , +4 , 5E-6+375 . ++2 ,5 .VA2 > 

OI< 6 >=1 0 . ++< 73 . 13-39 . 4+AL0GT< S >+5 . 5+ALOGK S >+ALOGT< 3 > > 

A 1=500, 

IF< 2-500 . > 49,60,60 
49 A 1 -Z 

6 0 FA 1 =9 , 8 06657< < 1 . + A 1 76 . 356766E+3 >++2 > 

FA 1 =FA 1 7< TX+ A2+ AT AH< T 1 +< A 1 - 1 25 . >+< 1 , +4 . 5E-6+< A 1 - 1 25 , >+*2 . 5 >/A2 > > 
D1=Z 

I F< 2-500/ 9 61,62,62 

61 D1 *500 . 

62 F0 1 =9 . 8 066 57 < < 1 , +0 1 76 . 356766E+3 >++2 > 

FD 1 =F0 1 /< TX+A2*ATAN< T 1 *< D 1 - 1 25 . >+< 1 . +4, 5E-6+< D1-1 25, >++2 .5 >/A2 > > 
N=0 

PREV=0 

SONE=< 0 1 -A 1 >*< FA t +FD1 >72. 

91 N=N+1 

IF <N-M> 92,92,95 

92 NINT as 2++N 
STWO=0. 

DEL=<01-At >7FLOAT< HINT > 

DO 93 1=1 , NINT, 2 
X 1 =A 1 +DEL+FLOAT < I > 

FX1 =9 . 80665/< < 1 . +X 1 76 , 356766E+3 >++2 > 

FX 1 =FX 1 7< TX+A2+ATA N< T 1 *< X 1 - 1 25 . >+< 1 . +4 . 5E-6+< X 1 - 1 25 . >++2 . 5 .VA2 > > 

93 STUIO-STWO+FXt 
CUR*=SQNE+4 , +DEL+STWO 

I F < EPS+ABS< CUR >-AB8< CUR-PREV > > 94 , 95 , 95 

94 PREV«CUR 

SONE»< SOHE+CUR >74, 

GO TO 91 

95 R=CUR/3 . 

TZ=TX+A2+ATAN< T 1 •+< Z-* 1 25 . >+< 1 , +4 . 5E-6+< Z- 1 25 . >++2 . 5 >7A2 > 

IF< 2-500. > 63,64,64 

63 R=-R 

64 OIT< 6>DI< 6 >+< S7T.Z >*EXP< -EI< 6 >+R7FK > 

C 

C++ LOOP BACK FOR ADDITIONAL CALCULATIONS 
C 

GO TO 56 
C 

C++ RETURN TO CALLING PROGRAM 

C 

999 RETURN 
END 

SUBROUTINE SLV< DEN, ALT, XLAT, DAY > 

C ++' +++++ + + +++ + + +++ + +++++++++++++++++++++*++++++++++++++*+++ 

C++ SUBROUTINE 'SLV' COMPUTES THE SEASONAL-LATITUDINAL ++ 



TABLE 1. (Continued) 


1019 C++ VARIATION OF DENSITY IN THE LOWER THERMOSPHERE < IN ** 

1020 C++ ACCORDANCE TO L, JACCHIA IN SAG 332, 1971. THIS AF- ++ 

1021 C++ FECTS THE DENSITIES BETWEEN 90 AND I60KM. THIS SUB- ++ 

1022 C ++ ROUTINE NEED NOT BE CALLED FOR DENSITIES ABOVE 160KM,++ 

1023 C++ BECAUSE NO EFFECT IS OBSERVED. +* 

1 024 C++ ++ 

1025 C++ THE VARIATION SHOULD BE COMPUTED AFTER THE CALCULA- ++ 

1026 C++ TION OF DENSITY DUE TO TEMPERATURE VARIATIONS AND THE++ 

1027 C++ DENSITY < DEN > MUST BE IN THE FORM OF A BASE 10 LOG » ++ 

1028 C++ HO ADJUSTMENTS ARE MADE TO THE TEMPERATURE OR CONSTI-++ 

1029 C++ TUENT NUMBER DENSITIES IN THE REGION AFFECTED BY THIS++ 


1030 

C++ 

VARIATION. 

** 

1031 

C++ 



** 

1 032 

C++ 

DEN 

* DENSITY < LOG1 0 > 

** 

1 033 

C++ 

ALT 

- ALTITUDE < KM > 

** 

1 034 

C++ 

XLAT 

* LATITUDE < RAD > 

** 

1 035 

C++ 

DAY 

■m DAY NUMBER 

** 

1 036 

C++ 



** 


1037 C +++ ** * ************ ************************* ****** ** ** ** * * * 
1 038 C 

1039 C++ INITIALISE DENSITY < DEN > » 0,0 

1 040 C 

1041 DEN *0.0 

1 042 C 

1043 C++ CHECK IF ALTITUDE EXCEEDS 1-60KM? 

1 044 C 

1045 IF <160 - ALT > 999,5,5 

1 046 C 

1047 C++ COMPUTE DENSITY CHANGE IN LOWER THERMOSPHERE 

1 048 C 

1049 5 2 * ALT - 90. 

1050 X « -0,0013+Z+Z 

1051 Y * 0 . 01 72+DAY +1.72 

1052 P * SIN<Y> 

1053 SP « SIN< XLAT > 

1 054 SP » SP+SP 

1055 S * 0, 014*Z*EXP< ><> 

1 056 D * S+P+SP 

1 057 C 

1058 C++ CHECK TO COMPUTE ABSOLUTE VALUE OF 'XLAT' 

1 059 C 

1 060 IF < XLAT > 10, 15, 15 

1061 10 D ** -D 

1062 15 DEN « D 

1 063 C 

1064 C++ RETURN TO CALLING PROGRAM 
1 065 C 

1 066 999 RETURN 

1 067 END 

1 068 SUBROUT I NE SL VH< DEN , DENHE , XL AT , SD A > 

1 069 C+++ **************************************** ************** 

1070 C++ SUBROUTINE 'SLVH' COMPUTES THE SEASONAL-LATITUDINAL ++ 


1071 C++ VARIATION OF THE HELIUM NUMBER DENSITY < ACCORDING ++ 

1072 C++ TO L. JACCHIA IN SAO 332, 1971 >. THIS CORRECTION ++ 

1073 C++ IS NOT IMPORTANT BELOW ABOUT 500KM. ++ 

1 074 C++ ++ 

1075 C++ DEN « DENSITY <LOG1 0> ++ 

1076 C++ DENHE * HELIUM NUMBER DENSITY < LOGIC) ++ 

1077 C++ XLAT « LATITUDE < RAD > ++ 

1078 C++ SDA « SOLAR DECLINATION ANGLE < RAD > ++ 


28 



TABLE 1. (Concluded) 


1 079 
1 080 
1 081 
1 032 
1 083 
1 084 
1 085 
1 086 
1 087 
1 088 
J 089 
1 090 
1 091 
1 092 
1 093 
1 094 
1 095 
1 096 
1 097 
1 098 
1 099 
1 1 00 
1 1 01 
1 1 02 
1103 
1 1 04 

I 1 05 

I I 06 
1 1 07 
1 1 08 

I 1 09 
1110 

II 1 1 
1 1 12 

1113 

1114 

1115 
1 1 16 

1117 

1118 

1119 

1120 
1121 
1 122 
1 123 
1 124 
1 1 25 
1 126 

1127 

1128 
1129 
1 130 
1131 


C>r::** ** ** **■•*•. * ** * ffc * * * * * * + * Ht *** sfc:*** sfc* sfc *:4< 

C 

C** PERFORM CALCULATIONS 
C 

DO » 1 0 , **DENHE 
A * 0,65*<SDA/0,40909079> 

C 

C** CHECK TO COMPUTE ABSOLUTE VALUE OF 'A ' 

C 

I F< A > 5,10,10 
FA * -A 
10 B » 0 , 5*XLAT 
C 

C** CHECK TO COMPUTE ABSOLUTE VALUE OF '& ' 

C 

I F< SD A > 15,20,20 
15 B - -B 

C 

C** COMPUTE X,Y,DHE, AND DENH E 
C 

20 X * 0 * 7854 - B 
V « SIN< X > 
y m y*y*y 
DHE» A*<Y ~ 0, 35356 > 

DENHE « DENHE + DHE 
C 

C** COMPUTE HELIUM NUMBER DENSITY CHANGE 
C 

D 1 « 10.**DENHE 

DEL= D1 - DO 

RHO* 1 0 , •♦.♦DEN 

DR HO - <6.646E-24>*DEL 

RHO * RHO + DRHO 

DEN * ALOGT< RHO > 

C 

C** RETURN TO CALLING PROGRAM 
C 

RETURN 

END 

BLOCK DATA 

C+* THIS IS THE BLOCK OATH FOR THE ** 

C** 'J7QMM ' PROGRAM. ** 

c*** ** ; 4« -b it* >♦< * :* ><< * H<?ft :*of: :* + A<>b * * Jfc :+t * 

C 

C*+ COMMON BLOCK DATA 
C 

COMMON XEDATAX IVR, IDA, MN, IHR, MIN, XM JO, F 1 0,F1 08, GI, ILAT, ILHG, 
XLAT, XLNG, I 1 

COMMON /FDAT AX A< 3 00 , 6 > , DENLG< 3 0 0 > , DENS< 3 0 0 > , DL< 3 0 0 > , EM< 3009, 

. T2Z< 300 >, XLATT< 300 >, XLONGC 300 >, 2< 300), XTEMP< 300 > 

END 
END* 
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Number of Data Values: 252 Mean Value? ,2522E 11 



TABLE 3. DESIGN SPACE STATION G&C SYSTEM MEAN 
TOTAL DENSITY FOR A LOW INCLINATION ORBIT 


Orbital Altitude* Total Density (kg/m 3 ) ** 


1100 

km 

(594 

n.mi.) 

0.5189 

X 

10 13 

1000 

km 

(540 

n.mi.) 

0.1018 

X 

1C' 12 

900 

km 

(486 

n.mi.) 

0.2105 

X 

10 12 

800 km 

(432 

n.mi.) 

0.4567 

X 

10~ 12 

700 

km 

(378 

n.mi.) 

0.1042 

X 

10~ U 

600 

km 

(324 n.mi.) 

0.2522 

X 

io" 11 

555 

km 

(300 

n.mi.) 

0.3847 

X 

10 " 11 

500 

km 

(270 

n.mi. ) 

0.6617 

X 

10" 11 

463 

km 

(250 

n.mi.) 

0.9719 

X 

10" 11 

400 

km 

(216 

n.mi.) 

0.1958 

X 

10" 10 


Ref: F 1q ? (230), A p (400), March 21, 1400 UT 


*The values of total density pertain to the 
kilometer altitude levels given. Nautical 
mile altitudes are approximated. 

**Design density values obtained from 
Reference 10. 



TABLE 4. 

PERCENTILES OF 

13-MONTH 

SMOOTHED 

VALUES OF 



F 

r 10. 7 

cm SOLAR RADIO 

FLUX AND 

A 

P 

GEOMAGNETIC 




INDEX 

OVER THE MEAN SOLAR 

CYCLE 






F10.7 




a p 


YR 

MO 

97.5% 

50% 

2.5% 


97.5% 

50% 

2.5% 

1992, 

05 

242.9 

154.0 

100.6 


18.1 

12.6 

10.5 


06 

239.0 

152.5 

96.2 


18.6 

12.7 

11.2 


07 

234.2 

150.5 

96.1 


19.1 

12.6 

11.3 


08 

230.8 

148.8 

95.6 


19.5 

12.9 

11.2 


09 

231.2 

143.3 

93.1 


20.1 

13.2 

11.4 


10 

230.7 

141.7 

90.9 


20.8 

13.3 

11.6 


11 

228.2 

140.2 

89.8 


21.2 

13.4 

11.7 


12 

226.2 

138.3 

92.2 


21.0 

13.9 

11.6 

1993, 

01 

225.7 

136.4 

89.7 


20.3 

13.8 

11.8 


02 

224.5 

134.5 

87.4 


20.5 

14.1 

12.0 


03 

223.2 

132.9 

89.9 


21.7 

14.3 

12.1 


04 

219.9 

126.9 

88.8 


22.5 

13.9 

12.0 


05 

215.0 

125.6 

84.9 


22.5 

14.2 

11.5 


06 

211.4 

124.5 

85.0 


22.1 

14.2 

11.2 


07 

205.9 

124.1 

83.2 


21.5 

13.9 

11.1 


08 

200.8 

123.5 

84.8 


21.3 

13.7 

11.3 


09 

195.4 

125.9 

82.8 


21.9 

13.3 

11.1 


10 

189.9 

124.9 

81.6 


22.9 

13.5 

11.4 


11 

185.5 

123.8 

79.7 


23.3 

13.5 

11.4 


12 

180.8 

122.7 

78.3 


23.2 

13.5 

11.4 

1994, 

01 

176.4 

121.5 

77.0 


23.1 

13.5 

11.0 


02 

176.7 

117.2 

76.6 


22.9 

13.1 

10.7 


03 

178.1 

118.2 

75.5 


21.6 

13.0 

10.9 


04 

177.6 

116.7 

76.0 


20.1 

13.8 

11.1 


05 

175.5 

115.4 

75.3 


19.5 

14.0 

11.6 


06 

171.4 

116.8 

73.1 


19.6 

13.7 

11.6 


07 

165.8 

114.7 

73.8 


19.9 

13.4 

11.4 


08 

161.9 

112.7 

70.8 


20.3 

13.5 

11.3 


09 

160.8 

113.4 

71.1 


20.5 

13.6 

11.4 


10 

159.9 

109.5 

71.7 


20.8 

13.6 

11.4 


11 

158.1 

108.1 

70.5 


20.9 

13.5 

11.2 


12 

155.6 

104.2 

70.0 


21.4 

14.2 

11.1 

1995, 

01 

151.4 

100.7 

70.7 


22.0 

14.0 

11.5 


02 

145.9 

99.7 

70.4 


21.8 

13.3 

11.5 


03 

139.8 

96.8 

68.8 


21.8 

13.2 

11.4 


04 

133.8 

97.8 

69.9 


22.0 

13.0 

11.4 


05 

128.6 

98.4 

72.5 


22.3 

12.3 

11.1 


06 

126.3 

98.7 

72.5 


22.6 

12.1 

11.3 


07 

124.7 

96.1 

71.6 


23.3 

12.3 

11.0 


08 

122.2 

93.6 

71.5 


24.0 

12.2 

11.2 


09 

119.9 

91.4 

71.6 


24.4 

11.7 

11.2 


10 

118.4 

90.7 

71.4 


24.7 

11.7 

11.2 


11 

117.5 

88.6 

70.5 


24.6 

11.7 

11.2 


12 

117.5 

87.5 

70.0 


24.2 

12.2 

11.2 
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TABLE 4 . (Continued) 


F 10.7 Ap 


YR MO 

97 . 5 % 

50 % 

2 . 5 % 

97 . 5 % 

50 % 

2 . 5 % 

1996 , 01 

116.3 

86.3 

70.1 

23.4 

12.2 

11.8 

02 

116.5 

86.3 

70.5 

22.6 

12.7 

11.4 

03 

117.5 

86.6 

69,8 

21.9 

12.8 

11.2 

04 

117.9 

87.2 

69.7 

21.6 

12.2 

11.1 

05 

118.0 

86.3 

68.2 

21.2 

11.9 

11.2 

06 

117.7 

85,4 

68.5 

20.9 

12.0 

11.3 

07 

116.8 

83.0 

67.6 

20.3 

12.1 

10.5 

08 

115.5 

82.1 

68.4 

19.6 

12.4 

9.8 

09 

113.5 

82.6 

67.9 

19.6 

12.5 

9.6 

10 

109.2 

81.8 

67.4 

19.6 

12.7 

9.1 

11 

103.5 

81.2 

67.3 

19.4 

12.9 

9.1 

12 

98.2 

79.9 

67.4 

19.0 

12.8 

8.8 

1997 , 01 

98.1 

79.9 

67.0 

18.5 

12.6 

8.7 

02 

98.8 

78.3 

67.1 

17.8 

12.6 

9.3 

03 

100.0 

77.4 

67.1 

16.9 

12.6 

10.1 

04 

100.4 

76.9 

67.3 

16.4 

12.5 

10.5 

05 

98.3 

76.6 

66.7 

16.6 

12.6 

10.1 

06 

95.2 

76.3 

66.9 

16.8 

12.7 

9.8 

07 

92.3 

76.0 

67.1 

17.0 

12.5 

9.4 

08 

91.0 

75.6 

67.2 

17.3 

12.4 

8.9 

09 

91.4 

75.1 

66.9 

17.6 

12.1 

9.1 

10 

91.6 

73.8 

66.8 

17.5 

12.0 

9.3 

11 

91.2 

73.2 

66.9 

17.3 

11.7 

9.3 

12 

90.8 

72.6 

67.1 

16.8 

11.6 

9.3 

1998 , 01 

90.2 

72.9 

66.7 

16.0 

11.5 

9.1 

02 

89.4 

71.8 

67.2 

14.6 

11.3 

9.1 

03 

88.5 

71.6 

67.2 

13.5 

11.0 

9.1 

04 

87.5 

71.3 

67.1 

13.6 

10.6 

8.9 

05 

86.3 

70.5 

67.2 

13.3 

10.4 

8.5 

06 

84.7 

70.4 

67.4 

12.8 

10.0 

8.1 

07 

82.6 

70.4 

67.8 

12.6 

9.5 

7.9 

08 

80.0 

70.1 

67.7 

12.3 

9.6 

7.6 

09 

77.6 

70.3 

68 . 0 

11.7 

9.6 

7.6 

10 

76.7 

70.8 

67. 8 

11.2 

9.6 

7.6 

11 

76.4 

71.7 

67.7 

11.0 

9.6 

7.7 

12 

77.0 

71.0 

67.6 

11.0 

9.7 

7.7 

1999 , 01 

78.0 

70.8 

67.8 

11.1 

9.7 

7.9 

02 

79.5 

70.4 

67.5 

11.1 

9.8 

8.1 

03 

81.6 

70.8 

67.7 

11.1 

9.9 

8.1 

04 

83.4 

70.9 

67.4 

11.4 

10.0 

8.1 

05 

85.8 

70.5 

67.7 

12.0 

10.2 

8.0 

06 

88.7 

70.5 

68.0 

12.9 

9.9 

7.8 

07 

91.7 

70 . 7 

67.8 

13.5 

10.2 

8.0 

08 

94.6 

70.9 

68.7 

14.1 

10.4 

8.1 

09 

99.1 

69.8 

68.6 

14.8 

10.5 

8.4 

10 

104.6 

71.5 

68.5 

15.6 

10.9 

8.3 

11 

110.0 

72.1 

68.7 

16.3 

10.6 

8.4 

12 

114.5 

72.9 

69.2 

17.3 

11.4 

8.4 
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TABLE 4. (Concluded) 


F10.7 A p 


YR MO 

97.5% 

50% 

2.5% 

97.5% 

50% 

2.5% 

2000, 01 

118.5 

73.9 

67.8 

18.2 

11,5 

8.4 

02 

122.0 

74.9 

68.5 

18.5 

11.4 

8.6 

03 

125.3 

76.1 

69,3 

18.6 

11.6 

8.8 

04 

130.2 

77.3 

70.1 

18.5 

11.8 

8.9 

05 

135.0 

81.4 

71.1 

18.2 

11.8 

9.3 

06 

140.5 

83.0 

72.0 

18.0 

11.8 

9.7 

07 

147.3 

81.8 

72. 9 

18.3 

11.7 

9.5 

08 

153.2 

83.6 

70.9 

18. 3 

11.9 

9.2 

09 

157.1 

85.6 

72.6 

17.5 

12.1 

9.0 

10 

159.7 

91.0 

74.6 

17. 0 

12.8 

9.5 

11 

161.4 

93.1 

73.6 

17.3 

12.8 

9.3 

12 

163.8 

95.2 

75.6 

17.3 

12.4 

9.5 

2001, 01 

166.8 

97.3 

74.1 

18.4 

12.5 

9.8 

02 

170.6 

99.0 

75.2 

19.7 

12.8 

10.1 

03 

174.7 

100.9 

72.7 

19.8 

12.9 

9.9 

04 

178.9 

99.5 

74.2 

19,8 

12.9 

10.0 

05 

184.5 

101.5 

75.1 

19.9 

13.0 

10.3 

06 

188.9 

103.6 

72.6 

20.2 

13.1 

10.2 

07 

190.2 

106.1 

75,5 

20.6 

12.7 

10.4 

08 

190.2 

108.8 

75.5 

20.7 

12.7 

10.7 

09 

189.4 

111.4 

75.9 

20.8 

12.6 

10.6 

10 

190.5 

113.4 

74.8 

21.0 

13.4 

10.7 

11 

193.6 

111.5 

75.8 

21,4 

13.5 

10.3 

12 

196.7 

116.6 

76.5 

21.9 

13.5 

10.8 

2002, 01 

198.4 

118.0 

77.8 

22.0 

12.9 

11.0 

02 

203.4 

119.2 

81.0 

20.8 

12.9 

10.8 

03 

209.8 

120.7 

84.2 

19.9 

12.7 

10.8 

04 

213.2 

122.4 

85.2 

19.6 

12.3 

10.7 

05 

216.0 

128.4 

86.7 

19.2 

12.0 

10.6 

06 

220.0 

130.8 

88.4 

19. 1 

12.1 

10.7 

07 

225.5 

133.1 

89.1 

18.9 

12.2 

11.1 

08 

228.3 

135.4 

91.1 

18.7 

12.4 

10.8 

09 

229.8 

137.9 

94.1 

18.5 

12.5 

10.7 

10 

231.9 

140.2 

96.5 

18.5 

12.4 

10.6 

11 

233. 9 

142.5 

99.0 

18.2 

12.5 

10.3 

12 

237.3 

149.1 

100.7 

18.1 

12.6 

10.5 

2003, 01 

241.1 

151.3 

97.4 

18.6 

12.7 

11.2 

02 

243.0 

153.0 

99.0 

19.1 

12.6 

11.3 

03 

242.1 

154.6 

97. 8 

19.5 

12.9 

11.2 

04 

241.2 

156.0 

100.7 

20.1 

13.2 

11.4 
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1 10 1 10 2 
USING MSFC— J70— 3 


10 3 10 4 10 5 10 6 10 7 10 8 10 9 10 10 10 11 10 12 1013 10 14 

CONSTITUENT CONCENTRATION (cm 3 ) 


FOR: MIN SOLAR CONDITIONS: 0400 hrs. USING Fio.7 = 70 & Ap = 0 

MAX SOLAR CONDITIONS: 1400 hrs. USING Fio.7 = 230 & Ap = 35 


Figure 1. MSFC/J70 constituent number density for 
low and high solar conditions. 




PERCENTILE 
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Figure 2. Percentiles of 13-month smoothed values of F10.7 cm solar 
^ radio flux over the mean solar cycle. 
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